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■ Abstract: We investigate for a large class of nonlinear wave equations, which allow 

for shock wave formations, how these solutions behave when they are "PT-symmetrically 
deformed. For real solutions we find that they are transformed into peaked solutions 
with a discontinuity in the first derivative instead. The systems we investigate include 
the PT-symmetrically deformed inviscid Burgers equation recently studied by Bender and 
Feinberg, for which we show that it does not develop any shocks, but peaks instead. In this 
case we exploit the rare fact that the "PT-deformation can be provided by an explicit map 
found by Curtright and Fairlie together with the property that the undeformed equation 
can be solved by the method of characteristics. We generalise the map and observe this 
type of behaviour for all integer values of the deformation parameter e. The peaks are 
formed as a result of mapping the multi-valued self-avoiding shock profile to a multi- 
valued self-crossing function by means of the 'PT-deformation. For some deformation 
parameters we also investigate the deformation of complex solutions and demonstrate 
that in this case the deformation mechanism leads to discontinuties. 
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1. Introduction 

Since the proposal that complex TT-symmetric quantum mechanical Hamiltonians may 
be viewed as self-consistent descriptions of physical systems 0], TT-symmetry has been 
exploited to propose and study many more complex extensions of real systems. Numerous 
new quantum mechanical and quantum field theoretical models have been investigated, see 
for instance ||, ||, || for recent reviews. Inspired by this success, the construction principle 
has also been used to suggest new classical models, such as complex extensions of stan- 
dard one particle real quantum mechanical potentials [||, |6|, [5], [8|, ||, [10| ], non-Hamiltonian 
dynamical systems [11], chaotic systems |l2j and deformations of many-particle systems 



such as Calogero-Moser-Sutherland models |13|, [14|, l5l 16, 17, 18, IS]. Here we will mainly 
focus on extensions of nonlinear wave type, such as the prototype Korteweg-deVries (KdV) 



equation [20, 21, |2^[ and closely related models [23, [25|, |26|, 27 1. In [27] it was demon- 
strated that when reduced these sytems can also shed light on large classes of one particle 
quantum mechanical models. 
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As is well known from transformation theory a lot of new information on nonlinear wave 
equations can be obtained from transformed equations by means of for instance Hopf-Cole, 
Miura or Backlund type. In a somewhat similar spirit we also exploit here the knowledge 
of a transformation in form of an explitly known PT-symmetrical deformation. Here we 
will mainly focus on the question of what kind of effect a PT-deformation has on a shock 
wave. It is well known that a shock forms when the crest of a wave overtakes the troughs. 
The challenge for a mathematical description is that one can no longer describe this phe- 
nomenon by a function since the surface of the wave becomes multi-valued. For real wave 
equations solvable with the method of characteristics, this happens when two character- 
istics cross each other or more generally when the first derivative becomes infinite. For 
PT-deformed equations remaining real, we argue here that the first derivative is discontin- 
uous but remains finite, whereas the second derivative tends to infinity. Solitonic solutions 



with this type of behaviour are often referred to as peakons [28|. When the deformation 
parameters are non odd integers, one is forced to consider complex solutions even in the 
undeformed case if one demands a real soltution for the deformed one. However, evolving 
this real deformed solution in time will convert it into a complex one. In addition, when 
compared to the real scenario, the peaks vanish and we observe discontinuties, which are 
generated due to the imposition of physical asymptotic boundary conditions. 

Our manuscript is organised as follows: In section 2 we describe the PT-deformed 
models we are investigating and how the explicit knowledge of the deformation map can 
be utilised to extract information about the systems, in particular the shock time and 
conservation laws. In section 3 we describe the general mechanism of how real shock waves 
are mapped into peaks and in section 4 how a modification of this mechanism leads to jumps 
in a complex scenario. In section 5 we present various numerical case studies supporting 
and illustrating our findings. We present our conclusions in section 6. 

2. Shock times and conservations laws from PT-deformation maps 

The model we wish to study here with regard to shock wave and peak formation is a 
PT-symmetrical deformation of a nonlinear wave equation 

PT £ : w t + f(w)w x = -> u t -if(u)(iu x ) £ = 0, (2.1) 

with f(w) being a well behaved function of w. Clearly the undeformed equation is PT- 
symmetric, that is being invariant under a simultaneous reflection in time t — > — t and 
space x — > —x, when PT : f{w) —> f(w), vu —> cw with c 6 C. The deformed equation is 



constructed as usual by taking into account that the PT-transformation is antilinear [29] 
and an overall minus sign in the second term is generated from i — > — i rather than from 
x — > —x. The special case f{w) = w corresponds to the PT-symmetrical deformation of 
the inviscid Burgers equation, also referred to as Riemann-Hopf or Euler-Monge equation, 
recently studied by Bender, Feinberg B3] and Curtright, Fairlie [S3]. The real version of 



the first equation in (2.1), especially for f(w) = w, appears mainly in fluid mechanics 
whereas its complex version is frequently encountered in the treatment of large N matrix 
models, see e.g. 
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Most wave equations are only directly related to their PT-symmetric deformations in 
the limit e — > 1 and an explicit transformation between the two is not known otherwise. 



For instance, this is the case for the deformations of the KdV-equation 21, 27]. In 
contrast, in some rare cases the deformation map VT e constitutes an explicit transforma- 
tion from the undeformed to the deformed system. As pointed out by Curtright and Fairlie 



|24] for the system fl2.1|) with f(w) = w, the first equation in (|2.1|) converts into the second 



under the map w i— > eu(iu x ) e ■ We generalise this here to equations with arbitrary f(w) 
to 

VT' E : w^ef{u)(iu x ) £ -\ (2.2) 
w t + ww x = -> ut — if(u)(iu x ) e = 0. 

This is seen as follows: Defining u = g(v) it is easy to verify that the deformed equation 
in (|2.1|) is converted into the deformed inviscid Burgers equation for v 

v t - iv{iv x f = 0, (2.3) 

provided the constraint 

v = f[g(v)][g'{v)r 1 (2.4) 

holds. Since we know already that w •->• ev(iv x ) £ maps the inviscid Burgers equation to 
the deformation fl2.3| ), we can derive from this the general map ( |2.2| ) when using (^4|) and 
v x = u x /g'(v). Note that the explicit form for the function g(v) is only required when we 



want to discuss equation (2^) in relation to the equation involving u. Essential is here 
only its existence, which follows from the fact that ( p.4|) is separable. 

Instead of VT' e , which transforms the inviscid Burgers equation, we can also construct 
VT £ , although in that case we have to be more specific about f(w). Taking for instance 



f(w) = w n the first equation in (2.1) converts into the second under the map 



VT e ■ w ^ yeu(iu x y- x . (2.5) 



Apart from the i, for n = 1 this reduces to the map found by Curtright and Fairlie [24|. 

These explicit maps ( |2.2| ) and ( |2.5| ) can now be exploited to investigate properties of 
the deformed equations. We start by considering in more detail how the solutions of the 
u>-equation in ( |2.lD are mapped into solutions of the deformed system. It is well known that 
the undeformed equation can be solved by the method of characteristics. The characteristic, 
i.e. the curve in the xt-plane at which w(x,t) = w(xo,0) ='■ wq(xo) is conserved, acquires 
in that case the form 

x = f(w )t + x . (2.6) 

The gradient catastrophe occurs when two of these characteristics cross or equivalently 
when w x tends to inifinity. Considering therefore 

, dx Q w' (x ) 
w x = w (x )—= (2.7) 

1 + 1 dx n 
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we can read off the earliest time, that is the breaking or shock time t™, and the correspond- 
ing position xf from ( [2.6D for which this happens 



tf = min ^-1/^^Tj > and xf = f [w (xf n )] tf + xf n . (2.8) 

Our concern here is how this translates into the deformed set of equations, i.e. what are 
the corresponding times £" and positions and moreover do the deformed systems exhibit 
shocks? 

is 



When f(w) = w n the shock time resulting from (2 

C = - ' 



(2.9) 



For n = 1 the expression Q2.9| ) agrees precisely with the formula (22) in [23], derived by 
analysing directly the deformed equation in ( |2.1| ) with the more complicated method of 
characteristic strips. Clearly we need to demand that the time is real, which is guaranteed 
when we replace no — > i a uo with uq 6 M. and a = (4m ± l)n/e, m G 7L. Thus for certain 
combinations of e and n we loose the possibility of shock wave generation for real solutions 
of the deformed equation. Nonetheless, in these cases we have a correspondence between 
a real undeformed solution and a complex deformed one. However, we will demonstrate in 
section 4 that if we do not insist in the undeformed solution to be real a shock formation 
is indeed possible, in constrast to the claims in 

Other quantities of interest which may be obtained from the explicit "PT-maps are 
conserved quantities. A conservation law for the first equation in ( |2.1|) is simply derived 
by multiplying it with nf(w) K and subsequent re- arrangement 

[/HI* + [/M K+1 L = o. (2.io) 

Therefore for any asymptotically vanishing function f(w) and constant k 6 1} the 

quantities 

4H = / f[w(x,t)] K dx, (2.11) 



are conserved in time. Correspondingly, we find for the deformed system the transformed 
charges 

/oo 
f[ef(u){iu x f^fdx. (2.12) 
-oo 

We will make use of these conserved quantities below. 



3. Peak formation mechanisms from real shock waves 

Having confirmed the result for the expression of the shock times by Bender and Feinberg 
in an alternative simpler manner, we diviate, however, from their interpretation of this 
result. Unlike Bender and Feinberg we conclude that these times correspond in general 
not to a gradient but rather to a curvature catastrophe, i.e. the first derivative stays finite 
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whereas the second tends to infinity. We reason as follows: For n = 1 it follows obviously 
from (HH) that 

w x = ie(iu x ) £ ~ 2 [u x + (e - l)uu xx ] . (3.1) 

This means that a possible shock in the "PT-symmetrically deformed inviscid Burgers 
equation, u x — > oo, would always correspond to a shock in the undeformed equation, that 
is w x — > oo. However, the reverse does not necessarily follow as a shock in the undeformed 
equation might correspond to u xx — > oo with finite u x , rather than to u x — > oo. Thus in 
the former scenario the shock time for the undeformed equation in ( [2.1D would correspond 
to a formation time of a different type of wave profile in the deformed equation in ([2.1|). 
We can identify the explicit form by expressing u in terms of w and use these expressions 
for our analysis. We find 



u (x,t) = (-iy--(s-i) 



£ e 



w(q, t) e- 1 dq 



w(q, t) s ~ 1 dq 



Differentiating this twice we obtain 

u x (x,t) = (— i) 1- e(e — l)^£~^w(x,t)~> 

u xx (x,t) = (-i) 1_ ^(e - l^^e^ (J w(q,t)^dqj ~w{x,l)- 
: j I — ;+ / w(q, i)~ dq ) w x (x, t)w(x, t)~ 



(3.2) 

(3.3) 
(3.4) 



which demonstrates that a shock in the undeformed ui-system will lead to u xx — > oo, as it 
directly depends on w x . On the other hand u x only depends on w. 

A closer inspection of the transformation (|3.2| ) explains how a shock is converted into 
a peak by means of the 'PT-deformation. Given the form of a shock profile for real w(x, t), 
as for instance depicted in figure [l] panel (a), we first need to convert the multi- valued 
profile into a single valued function, which is achieved by parameterising the profile by the 
arc length s. From (|3.2D we obtain 



u(s, t) := iu(x, t)~' 



1) 



, ._L_ dq 
» M -pds, 
ds 



(3.5) 



with arc length element ds = Vdw 2 + dx 2 . We compute u(s, t) for e = 3, where we take 
the positive square root for s < S3 and the negative one for s > S3. When transforming 
back from s — > x this function becomes multi-valued crossing itself as depicted in figure || 
panel (b) for a time t\ > t s . 

This selection of the branches will produce a peaked function. In general the choice 
of the different branches is naturally governed by the appropriate boundary conditions 
matching the initial profile. Note that we can eliminate the part s± — > S2 — > S3 — > S4 without 
destroying the consistency of the model and thus convert a solution from a multivalued 
one into a single valued peaked function. It follows by ( [2.11 ) that the conserved quantity 
I K with k = (e — 1) _1 is unaffected by this change and remains preserved, since for that 
choice f 00 = f Sl + f . However, for different values of n the I K no longer constitute 
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charges for the peaked solution. This argument is similar to the standard introduction of 
a shock front, the position of which is usually chosen in such a way that I\ is preserved. In 
principle this could be implemented for the undeformed system. It is clear that changing 
from u(x,t) to u(s,t) will not alter the argument very much, apart from introducing yet 
more possible branches. The peaked solutions are to be understood in the weak sense, such 
that a rigorous treatment requires the use of test functions. 




Figure 1: (a) A multi- valued self-avoiding shock wave w(x,ti). (b) Transformed multi- valued 
self-crossing shock wave u(x, ti) for s = 3. 

We have argued above that peaked solutions are the most common ones appearing in 
the deformed equations. However, under some special circumstances, that means particular 
initial profiles, we may also encounter shocks in the deformed system. We observe from 
Q3.3| ) that u x — > oo whenever J x w(q, dq = for e > 0. In turn this means by ( |3.2| ) 
that u — > 0. Indeed one may construct such type of solutions, as we will demonstrate 
below. 



4. Jump formation from complex shocks 

So far we have focused on real solutions w(x, t) for the inviscid Burgers equation. However, 
if we wish to consider a real transformed solution u(x,t) in the case for e not being an odd 
integer then in the light of (^T]) we are forced to consider also complex initial conditions. 
In fact, the complex inviscid Burgers equation has appeared before in the literature in 
several different contexts, see for example [31], [32] or p3] for applications to geostrophic 



flows or large N matrix models, respectively. 

The method of characteristics can readily be adapted to this case. Exactly like in 
the real scenario, the solution is given by w(x,t) = wo(xo(x,t)), where xo(x,t) is found 
by inverting fl2.6| ). Notice that in general xo(x,t) £ C even for real (x,t), so that the 
solution is defined by analytically continuing the initial condition to complex values of the 
argument. 



In [34 1 it has been revealed that the shocks in inviscid Burgers equation are due to the 



presence of square root singularities of w(x, t) in the complex x plane. The position of 
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the singularities at time t can be found by first determining xo by means of ( |2.7| ), and 
subsequently computing x = xq — /(^o(^o))/ ■ 

For t = all such singularities are located at complex infinity. Considering the solutions 
as functions of real x and t, with t being the parameter governing the flow, the singularities 
move in the complex x-plane, exhibiting a shock whenever they reach the real axis. From 
this discussion follows that in order to find all the possible shock times we need to impose 
two conditions on xq: 



Im 



df(w (x )) 

dxn 



and Im 



xo - f(w (x ))/ 



df(w Q (x Q )) 

dxn 



0. 



(4.1) 



There is a crucial difference between the real and the complex case. When wq, /(ioo) £ 
R, every xq G R solves the second equation in ( |4.1| ) and we encounter a shock for all times 
in some interval [t s , oo). The reason for this is that in the real case the singularities reach 
the real axis in complex conjugate pairs and thereafter, i.e. for t > t s , never leave it. In 
contrast, in the complex case the solutions of the second equation in ( [4.1[ ) are in general 
isolated points in the complex plane, so that a gradient catastrophe will be an isolated event 
in time. The mechanism responsible for the jump is the application of a matching condition 
between the intial boundary condition and the ones for the evolved solultion. Usually the 
boundary conditions are taken to be physical, that is asymptotically vanishing. Unlike as 
in the real case, one can no longer stitch the two asymptotic solutions W1/2 —> ±oo together 
in a continuous manner, but instead one is forced to introduce a jump. 

We shall now support and illustrate our general findings with some numerical studies. 



5. Numerical case studies 



5.1 Real w and real u 

The e = 3 deformation with f(w) = w is the simplest example allowing for real solutions 
for the undeformed as well as for the deformed equation. In order to be able to compare 
directly with the results in [^3| we consider here the same initial profile of the form of a 
Cauchy distribution uq = (1 + X 2 )" 1 , such that by ( ^ ) the corresponding initial profile in 
the undeformed equation results to wq = — 12x 2 /(l + x 2 ) 5 . According to (|2.9| ), the gradient 
catastrophe occurs therefore when 



C(*o) 



(i + *3 



2\6 



24x (4xjj - 1)' 

This function has two distinct minima, which we identify both as shock/peak times 



(5.1) 



1 



x, 



0,1 



X, 



0,2 



6V2 

-1 

6V2 



23 



/385 « 0.216621, = 
23 + \/385 ps -0.769392, tf 2 



(95 - v / 385) ( 



2 21 3 10 (5711-2^/35; 
(95 + \/385) 6 (5 + ^385 



0.311791, 



(5.2) 



0.644466. 



2 18 3 10 a/2(23 + y / 385) 
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In a system with shocks the second time is usually not easy to realise numerically, but 
for the deformed systems these times correspond to peaks and are directly accessible in 
numerical simulations. The numerical values for the times agree with those provided in 
23]. In addition we compute from ( p.§| ) the corresponding positions of the shocks/peaks 



as 



365) 



64 5V 11 - 2V35 



0.0770263, 



23 + V385 



L s,2 



16 



-1.21712. (5.3) 



Figure || panel (a) exhibits how a shock develops in the undeformed system at the time 
tf 1 and position xf^ as predicted by (|5.2|) and (|5.3|), respectively. In panel (b) we observe 
that a peak develops at the same times and positions = t™ 1} x"^ = xf 1 and also 
at tp 2 = tf 2 i x s2 = x sV For the- deformed system the numerical integration over the 
discontinuities does not pose any major obstacle and the event of the second peak can be 
simulated simply by integrating until that time. 




Figure 2: (a) Solution of the inviscid Burgers equation for transformed Cauchy distribution initial 
profile, (b) Solutions of the e = 3-deformation with Cauchy distribution initial profile. 

In contrast, for the undeformed system this is not possible in such a straightforward 
manner because after the first shock time the function w(x,t) becomes multi-valued, such 
that the standard procedure becomes meaningless. In principle this problem can be over- 
come by introducing a shock front and the preservation of some conserved quantities as 
argued above. For a detailed survey on these techniques see e.g. |35|. Here this is not 
needed and instead we can use the fact that our undeformed system is implicitly solved 
by w = u>o(x — wt) which we can solve numerically for w(x,t). Subsequently u(x,t) is 
computed from (p^). 

The results of this computation are depicted in figure |3], which shows that the second 
shock time i Sj 2 and position x s ^ are approached at their predicted values ( |5.2[ ) and (|5.3p, 
respectively. These positions in space and time coincide with those of the second peak in 
the deformed system. In panel (b) we observe that the effect on the smoothness of the 
curve is much less pronounced for the second "peak". In this case it is hardly visible due 
to fact that it is not located on the crest of the wave. 
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Figure 3: (a) Solution of the inviscid Burgers equation at times before and after the second shock 
formation with t = 0.4 dotted (red), the second shock time t = t s ^ solid (black) and t = 1.0 dashed 
(blue) for transformed Cauchy distribution initial profile, (b) Solution of the e = 3-deformation of 
the inviscid Burgers equation at the second shock time t = t s ^ for transformed Cauchy distribution 
initial profile. 



Let us now present an example for the formation of a shock rtaher than a peak. For 
the deformed system we take as initial profile uq = x/(l + x 2 ), such that the initial profile 
of undeformed equation results to wo = — 3x(l — x 2 ) 2 /{l + x 2 ) 5 . According to (2JS) this 
system develops a shock at time t™ = 1/3 at position xf = 0. For the undeformed system 
these findings are clearly confirmed by our numerical results depicted in figure ||] panel (a). 
Notice in panel (b) that when the solution evolves beyond the shock time two more shocks, 
in the sense that u x — > oo, develop when u becomes zero. This behaviour is forced by the 
relation between w and u, (EO), since w is finite. 





Figure 4: (a) Shock wave formation at w — for the inviscid Burgers equation at times t = 
dotted (red), the shock time t s = 1/3 solid (black), t — 1 dashed (blue) and t = 2 dasheddotted 
(green) for the transformed initial profile x / (1 + x 2 ) . (b) Solution of the e = 3-dcformation of the 
inviscid Burgers equation at u = 0. 
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5.2 Real w and complex u 

The e = 2 deformation with f(w) = w is the easiest case to investigate the scenario for 
which w £ R and u ^ R. We compute initially the time and position of the shock and 



peaks in the undeformed system (2.1), respectively. 




Figure 5: (a) Solutions of the inviscid Burgers equation for a transformed Gaussian initial profile, 
(b) Solution of the e = 2-deformation with a Gaussian initial profile in the exp(i7r/4)-direction. 



We take the initial profile to be of Gaussian form uq 



! 4 leading to wq 



— Axe~ with the help of the transformation (|2.5|) . As argued above, to take the phase in 
uq is one possibility to guarantee real shock times. The time for the gradient catastrophe 
computed from (2.9) is t™ c = e~ x o/(A — 16x 2 ), which becomes minimal for x™ 111 = 0. The 
resulting shock/peak time and position are therefore tf = = 1/4 and xf = = 0, 
respectively. 

Figure [5] panel (a) clearly shows how a shock develops in the undeformed system at the 
predicted shock time t s and position x s . Panel (b) exhibits that in the deformed system this 
shock is converted into a peak occurring at the same time and position. In general this wave 
travels in the complex u-plane, but here we have only plotted here the exp(i-7r/4)-direction 
for which u(x,t) becomes real. We also note that the peak becomes more pronounced as 
the wave evolves beyond its time of formation. 



5.3 Complex w and complex u 

Here we will consider the deformation with e = 3/2 . As discussed in section |4], we now 
require a complex solution for the undeformed equation in order to generate shock waves. 
We select a simple initial condition which vanishes asymptotically 



w (x) 



w(x, 0) 



e 4 



x 2 + 1 



(5.4) 



and in addition allows to solve the well known implicit realtion w = wq{x — wt). The 
transformation ( |2.2[ ) then guarantees that we have real initial conditions for the deformed 
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equation 

with boundary conditions 
lim Uq(x) = 0, 



Uq[X) 



lim no (a;) 

-i-+oo 



(y 2 + 1) 2 



dy 



1/3 



=: fc « 1.2794, and 



lim d x u Q (x) 

x— >±oo 



0. 



(5.5) 



(5.6) 



We demand that these boundary conditions are preserved for t > 0. 

Given wq, we compute the shock time as outlined in section Q The conditions ( |4.1| ) 
have two solutions: zo,i ~ 0.164903 — 0.553299i and zo,2 = — £o,l) an d> correspondingly, 
we find two shock times: t Sj i ~ 0.4791 > and t Sj 2 = —t s ,i < 0. The shock occurring for 
positive time is located at x Sj \ « 0.494709. 

In figure |6|, we plot real and imaginary part of the solution w(x,t), obtained by inverting 
the relation w = wo(x — wt). In this example, w(x,t) has three branches, of which only two 
are represented in the figure. The branch represented as a solid (black) line is the solution 
satisfying the initial condition ( |5.4j ). 

We see that the two branches touch at t = t s i- For t > t s ±, the left part of one branch 
has connected with the right part of the other, imposing a jump on the physical solution 
in order to preserve the asymptotic behaviour lim x ^ l .± 00 w(x,t) = 0. 

To see how this reflects on the evolution of the deformed field, we have constructed 



the solution u(x,t) using the relation (3.2) 



u(x, t) 



u r x . , 



dy 



1/3 



(5.7) 



and evaluated the integral numerically. In figure fj], we represent u(x,t) for a sequence of 
times leading to t Sj i, and we observe again that u x is continuous, while u xx {x Sj \) — > oo as 
t — > t s i. 



After the shock time, w(x, t) develops a jump. We will show that the boundary condi- 
tions (|5.6| ) impose a jump on u(x,t) as well. In order to see this we will try to match the 
following expressions: 

u(x,t)= (~ j X w^(y,t) 2 dyys; (5.8) 

u(x, t) =(k 3 -j f + (y,t) 2 dy) (5.9) 

where the two branches tf/ 1 ' 2 ) (x, t) are defined by: w^(— oo,t) = = (+oo, t) for 
t > t S} i- By construction, u(x,t) satisfies the left boundary condition in fl5.6j ), while u(x,t) 
satisfies the condition on the right. 

In figure |8] we show u(x,t) and u(x,t) for t = 1. We notice that, contrary to the real 
case treated in section 3, we can not find an x* such that n(x*,t) = u(x*,i). This can 
be understood because the former condition now splits into two real equations, while we 
have only one real parameter to tune. The consequence is that a continuous solution 
for u(x,t) does not exist beyond t s> i. 
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Figure 6: (a), (b) Real and imaginary part of the solution of the inviscid Burgers equation with 
complex initial condition (5.4) solid (black) and a second, unphysical branch dashed (red) at time 
t = 0.45 < t a i. (c), (d) The two branches touching at the shock time t = t St \. (e), (f) The 
physical solution with vanishing asymptotic boundary conditions exhibits a jump discontinuity at 
t = 0.55 > t s i. 



5.4 Deformations with odd g and f(w) = w 

We will now study how the systems behave as a function of increasing values of the defor- 
mation parameter e. We keep the initial profile in the deformed system to be a Cauchy 
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Figure 8: Real and imaginary parts of the two functions u(x, t) solid (black) and u(x, t) dotted (red) 

1. The condition Re u(x,t) = Re u(x,t) is satisfied for x — x\ ~ 1.0663, 
Im u{x 1 1) for x — X2 ~ 0.1893 ^ x\. 



defined in (5.8), for t 
and the condition Im u(x,t) 



distribution uo = (l+x 2 )^ 1 , such that by (p.5| ) the initial profile in the undeformed equa- 
tion will change. We do not report the explicit expressions for these functions here, but 
only the resulting shock and peak times the following table: 



e 




t s ,2 


x s ,i 


X s ,2 


3 


0.311791 


0.644466 


0.0770262 


-1.21712 


5 


0.394011 


0.662872 


-0.18255 


1.05226 


7 


0.594697 


0.913866 


0.241058 


-0.970114 


9 


0.997223 


1.45053 


-0.279227 


0.919109 


11 


1.78617 


2.50127 


0.306641 


-0.883621 


13 


3.34619 


4.555 


-0.327569 


0.857142 



In figure ^ panels (a), (b) we observe that the first shock times are reproduced in our 
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numerical solutions and that they occur at the predicted locations. 




-2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 

X X 



Figure 9: (a), (b) Solutions of the undeformed system at the shock time i s ,i for various values of 
e with transformed Cauchy distribution initial profile, (c), (d) Solution of the deformed system at 
the shock time t Sj i for a Cauchy distribution initial profile. 

We observe that the overall qualitative behaviour does not change with increasing values 
of e, but the amplitutes are considerably reduced. We also note that the systems with 
e = 4m — 1 and e = 4m + 1 with m G N break first on the right and left wave front, 
respectively. In the panels (c) and (d) we observe that the shocks have been smoothed out 
considerably and that the larger e becomes the less pronounced they are. The peaks are 
extremely small when they occur on the tails of the waves rather than on its crest. 

5.5 Multi-peak solutions 

In the same way as for as for peakons an interesting question is whether it is possible to 
obtain multi-peaked solutions. We demonstrate here that the answer to this is affirmative. 
As an example we consider the equations with f(w) = w and deformation parameter e = 3. 
Taking then for instance the initial profile in the deformed system to be a sum of two shifted 
Cauchy distributions no = [1 + (x— l) 2 ] -1 + [1 + (x+ 1) 2 ] -1 , the initial profile of undeformed 
equation results to wq = —96 [x 2 + 2) (x 5 + 4x 3 — 4x) 2 (x 4 + 4) 5 . By means of ( |2.9D the 
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time for the gradient catastrophe is then computed to 

(x 4 + 4) 6 

t w (xn) = I 1—1 f5 10) 

gcV ; 384x (2x u + 25x 12 + 60x w - 156x 8 - 384x 6 + 240x 4 + 192a; 2 - 64) ' v ' ; 

Minimising this function we find four shock/peak times t s \ = 0.221045, t S 2 = 0.429609, 
t s3 = 0.558845, t s4 = 0.798264 and corresponding positions x sl = 1.01299, x s2 = -2.21359, 



x S 3 = —0.856069, x S 4 = 0.116185. All these values are accurately reproduced in figure |10 
Once again the peaks emerging on the crest of the wave are well pronounced. 




IE 



Figure 10: (a) Solutions of the inviscid Burgers equation for transformed sum of shifted Cauchy 
distribution initial profile at shock peak times t s \ dotted (red), t s2 dashed (blue), t s3 dasheddotted 
(green) and t S 4 solid (black), (b) Solutions for its e = 3-deformation at the same times and Cauchy 
distribution initial profile. 



6. Conclusions 

For large a class of nonlinear wave equations of inviscid Burgers type ( p.l[ ) we have shown 
that their real shock wave solutions are smoothed out in their "PX-symmetrically deformed 
counterparts into peaks. The mechanism for the peak formation was identified to be the 
folding of a self-avoiding multi-valued function into self-crossing multi-valued function. 
Under the preservation of certain conservation laws one can consistently eliminate the 
looping part of the self-crossing multi-valued function and thus converting it into a single 
valued weak solution of the deformed wave equation. In general, we found that the larger 
the deformation parameter e the smoother the peaks. Our analytical arguments were 
facilitated by the explicit knowledge of the 'PT-transformation from one system to the 
other. 

We also showed that shocks in the complex solutions for the undeformed equations 
will lead to discontinuities in the solutions for the deformed equations. 

It would be of great interest to construct more explicit maps for different types of 
nonlinear wave equations to their PT-symmetrically deformed counterparts, such as for 
instance for the KdV equation. 
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As a practical application one may potentially use these observation in numerical inves- 
tions of nonlinear wave equations with shock wave formations. In general those equations 
are difficult to handle, even numerically, but one can exploit our obervations by first solving 
the deformed equations, which are simpler to deal with as they just exhibit peaks and then 
transform those solutions to the undeformed system. 

Acknowledgments: AC is supported by a City University Research Fellowship. 
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